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We study the relation between stochastic and continuous transport-limited growth models, which 
generalize conformal-mapping formulations of diffusion-limited aggregation (DLA) and viscous fin- 
gering, respectively. We derive a nonlinear integro-differential equation for the asymptotic shape 
(average conformal map) of stochastic aggregates, whose mean-field approximation is the corre- 
sponding continuous equation, where the interface moves at its local expected velocity. Our equation 
accurately describes advection-diffusion- limited aggregation (ADLA), and, due to nonlinear averag- 
ing over fluctuations, the average ADLA cluster is similar, but not identical, to an exact solution 
of the mean-field dynamics. Similar results should apply to all models in our class, thus explaining 
the known discrepancies between average DLA clusters and viscous fingers in a channel geometry. 
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Developing effective mean field approximations to non- 
linear stochastic equations constitutes a major challenge 
in various active fields of statistical physics, e.g. hydro- 
dynamic turbulence [1] and self organized criticality [2] . 
Straightforward derivation of such approximate theories 
typically involves uncontrolled assumptions required to 
"close" an infinite hierarchy of equations for moments 
of the probability distribution of the stochastic field. An 
alternative approach consists of deriving asymptotic solu- 
tions to a deterministic version of the original stochastic 
dynamics, assuming that such solutions capture the be- 
havior of ensemble average of the original stochastic field 
[2]. Such approach, however, might turn out to be unre- 
liable as well, since it underestimates the possible effects 
of noise on the asymptotic evolution of a stochastic field. 

A nontrivial example in which such approach has been 
advanced over the last two decades is the fractal mor- 
phology of patterns observed in computer simulations of 
the celebrated diffusion limited aggregation (DLA) model 
[3]. Since the relation between the mathematical formu- 
lations of the stochastic DLA process and the contin- 
uous viscous fingering dynamics was established [4], the 
striking similarity between patterns observed in both pro- 
cesses has triggered various attempts to interpret viscous 
fingering dynamics as a mean field for DLA [5-8] . 

In this Letter, we study the connection between a 
broad class of stochastic transport-limited aggregation 
processes and their continuous counterparts [9]. In our 
models, growth is fuelled by nonlinear, non-Laplacian 
transport processes, such as advection-diffusion and elec- 
trochemical conduction, which satisfy conformally invari- 
ant equations [10]. Stochastic and continuous dynamics 
are defined by generalizing conformal-mapping formula- 
tions of DLA [11] and viscous fingering [12, 13], respec- 
tively. We show that the continuous dynamics is a self- 
consistent mean-field approximation of the stochastic dy- 
namics, which, nevertheless, does not accurately predict 
the average shape of a random ensemble of aggregates. 

We consider a set of two-dimensional scalar fields, ip = 
{<px, tf2, ■ ■ ■ , <Pm}, whose gradients produce quasi-static, 
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FIG. 1: (a) A typical ADLA cluster of n = 10 5 "particles" (it- 
erated conformal maps). Color contours indicate the particle 
concentration c, and solid yellow curves are fluid streamlines, 
(b) The exact, asymptotic shape of the analogous continuous 
dynamics, G c (w), which describes solidification in a fluid flow. 
For both (a) and (b), we set Pe ~ 20 so that the far fields are 
similar. 

conserved "flux densities" , 

M 

Pi=S^( V )V¥>i, V-F s: = (1) 

in Q z (t), the exterior of a singly connected domain that 
represents a growing aggregate at time t. (The coeffi- 
cients, {Cij} may be nonlinear functions of the fields.) 
The crucial property of the nonlinear system (1) is its 
conformal invariance [10]: If ip(w,W), not necessarily har- 
monic, is a solution in a domain Q w and w = f(z) is a 
conformal map from Q z to O^, then y>(/(z), f(z)) is a so- 
lution in Q z . Using this fact, the evolving domain, il z (t), 
can be described by the conformal map, z = g(w, t), from 
the exterior of (say) the unit disk, Q w , 

Growth is driven by a combination of flux densi- 
ties, Q = ~^2^ = i Bi((f)Fi, on the boundary with a lo- 
cal growth rate, a = n ■ Q, where n is the unit nor- 
mal vector at z £ dQ. z (t). For continuous, deterministic 
growth, each boundary point z moves with a velocity, 
v(z) = acr(z) n{z), where a is a constant. For discrete, 
stochastic growth, the initial seed, fi z (io = 0), is itera- 
tively advanced by elementary "bump" maps represent- 
ing particles of area Aq at times t\, ■ ■ ■ , t n . The waiting 
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time t n — t n -i is an exponential random variable with a 
mean set by the total integrated flux. The probability 
density to add the nth particle in a boundary element 
(z,z + dz) e dQ z (t n -i) is proportional to a(z)\dz\. 

The classical models are recovered in the simplest case 
of one field (M = 1). DLA corresponds to stochastic 
growth by diffusion, F = Q = -DVc, from a distant 
source (c <~ log|z| as \z\ — > oo) to an absorbing cluster 
(c = for z G Cl z (t)), where c is the particle concentra- 
tion and D the diffusivity. Viscous fingering corresponds 
to continuous growth by the same process, where c be- 
comes the fluid pressure and D the permeability. 

The simplest, nontrivial models with multiple fields 
(M = 2) involve diffusion in a fluid flow. The 
stochastic case is advection-diffusion-limited aggregation 
(ADLA) [9], illustrated in Fig. la. Particles are de- 
posited around a circular seed of radius, L a , from poten- 
tial flow, v = V<^, of uniform velocity U far from the 
aggregate. The dimensionless transport problem is 

Pe V<^ • Vc = V 2 c, VV = 0, zefl z (t) (2) 
c = 0, n-Vip = 0, a = n-Vc, z e dQ. z {t) (3) 
c — ► 1, V<£ — ► x, \z\ — > cxo, (4) 

where c is the concentration of particles. Here x, ip, c, 
and a are in units of L Q , UL a , C, and DC/L , respec- 
tively, and Pe Q = UL Q /D is the initial Peclet number. 
Numerical solutions and asymptotic approximations are 
discussed in Ref. [14]. 

The transport problem is conformally invariant, except 
for the boundary condition, Eq. (4), which alters the flow 
speed upon conformal mapping. Instead, we choose to 
fix the mapped background flow and replace Pe Q with 
the renormalized Peclet number, Pe(t) = Ai(i)Pe , when 
Eq. (2) is transformed from £l z (t) to Q w . The "confor- 
mal radius", Ai(t), is the coefficient of w in the Laurent 
expansion of g(w,t) and scales with the radius of the 
growing object [11, 15]. Since A\{t) — ► oo for any initial 
condition, the flux approaches a self-similar form, 

cj(0;Pe) ~ 2 y / Pc/^ sin(0/2) as Pc oo. (5) 

More generally, there is a universal crossover from DLA 
(a = constant) to this stable fixed point, where Pe(t) = 
j4i(i)Pe is the appropriate scaling variable [16]. 

The continuous analog of ADLA is a simple model for 
solidification from a flowing melt [17]. More generally, 
continuous dynamics in our class is described by a non- 
linear equation, 

Re (w g' g t ) = a a(w; Pe(£)) for |io| = 1. (6) 

which generalizes the Polubarinova-Galin equation for 
Laplacian growth [12, 13] (a = 1). In the case of 
advection-diffusion [17], only low-Pe approximations are 
known, but we have found an exact high-Pe solution of 
the form, g(w,t) = A\(t) G c (w), where 

A 1 (t)=t 2 / 3 , G c (w)=wy/l-l/w. (7) 



This similarity solution to Eq. (6) with aa(6, t) — 
y/ Ai(t) sin(0/2) describes the long-time limit, according 
to Eq. (5). (We do not know the uniqueness or stability 
of this solution or whether it can be approached without 
singularities from general initial conditions.) Just as the 
Saffman- Taylor finger solution (for a — 1) has been com- 
pared to DLA in a channel geometry [18], this analytical 
result begs comparison with ADLA. 

As in Ref. [9], we grow ADLA clusters by a modi- 
fied Hastings-Levitov algorithm [11]. The random at- 
tachment of the nth particle to the cluster is described 
by perturbing the boundary dfl z (t n -i) — ► dQ z (t n ) by 
a "bump" of characteristic area Ao- This leads to the 
recursive dynamics 

9n(w) = g n -i o </>\ n ,e n (w), g n {w) = g(w,t n ) (8) 

where 4>\,e(w) is a specific map, conformal in Cl w , 
that slightly distorts the unit circle by a bump of 
area A around the angle 9. The parameter, A„ = 
^o\g' n -i( et0n )\~ 2 is the area of the pre-image of such 
bump under the inverse map g^ 1 . The angle 6 n is chosen 
with a probability density p(6;Pe(t n )) oc a(e l6 ; Pe(i n )), 
so the expected growth rate is the same as in the contin- 
uous dynamics. 

For an ensemble of n-particle aggregates, a natural def- 
inition of average cluster shape is the conformal map, 
(G n (w)), defined by averaging at a point w £ £l w all the 
maps, G n (w) = g n (w)/ Ai(t n ), rescaled to have a unit 
conformal radius. We then ask: What is the limiting 
average cluster shape, (G OCl (w)) — lim„^ 00 (G„(w)) ; and 
how does it compare to the similarity solution, G c (w), of 
the continuous growth equation (6)? The same questions 
apply to any of our transport-limited growth models, but 
here we focus on ADLA as a representative case. 

To provide numerical evidence, we grow 2000 ADLA 
clusters of size n = 10 5 using the semi-circular bump 
function in Ref. [15] (with a = 1/2). To reduce fluc- 
tuations, we aggregate small particles, Ao = 1/16, on 
a large initial seed (go(w) — w, \w\ = 1). To reach at 
the asymptotic limit faster and also match the assump- 
tion of G c (w), we fix the angular probability measure, 
Poo(9) = sin(0/2)/4 for Pc = oo, throughout the growth. 
In Fig. 2a, we plot the average contour of the ensemble, 
(G n {e ie )) at n = 10 5 along with that of the continuous 
solution, G c {e ie ). To give a sense of fluctuations, we also 
plot a "cloud" of points, G n {e 10 ) over uniformly sampled 
values of 9. Fig. 2b is the zoom-in of the boxed region 
in Fig. 2a, where we also show (G n (e ie )) at n = 10 3 
and n = 10 4 . Although the convergence of (G n (w)) is 
easily extrapolated, the n = 10 5 line has not reached 
at the asymptotic limit yet. The branch point at w = 1 
seems to be related to the slow convergence. Ignoring the 
unconverged area, (G n (w)) and G c (w) are quite similar, 
and yet clearly not same. The average, (G n (w)), better 
captures the ensemble morphology reflected by the cloud 
pattern than G c (w) and the opening angles at the branch 
point of the two curves are also different [16]. 
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FIG. 2: (a) The "cloud" of the sampled points, G n (e ie ) and its 
mean field contour, (G„{e lB )) at n — 10 5 (dashed) displayed 
with the asymptotic limit, (Goo{e' e )) and the steady state 
shape of the continuous growth, G c (e l8 ). (b) The zoom-in of 
the box in (a). Two contours, {G n (e l9 )}, at n = 10 3 (dotted) 
and 10 4 (dashed dot) are added to show the slow convergence 
to (Goo(e* )) near the rear stagnation point. 



Next we derive an equation for {G n (w)) in the asymp- 
totic regime. For growing aggregates A„ — ► as n — > 
oo [15]. Following Hastings [19], we use Eq. (8) to 
derive a linearized recursive equation for G n +\{w) for 
|to — e l8n+1 \ yJXn+x- Letting (A, 6) denote the param- 
eters of the (n + l)th bump, we obtain: 



G„+iH - (1 - a\)(G n (w) + a\H e (w)G' n (w)) 
~ G n (w) + a\(He(w)G'Jw) - G n (w)). 



(9) 



where Hg(w) — w (w + e l9 ) / (w — e l9 ) and we use 
Ai(t n+ i) = (1 + A) a Ax{t n ). 

Stationarity of the ensemble of rescaled clusters im- 
plies: 



(GnH) 



d9 Poo (0)(G n+1 {w)) . 



(10) 



Our analysis applies for conformally invariant transport- 
limited growth from an isolated seed with general angular 
probability distributions, although we will focus on the 
case of ADLA, Poa (9) = sin(0/2)/4. 

Using Eq. (9), we get the fixed-point condition: 



d9 Po o(9)(XG 00 (w)) 



d6 Poo (0)(XGM)H 8 (w). 

To facilitate further analysis, we approximate the left 
hand side of Eq. (11) as 



2tt 



d0 Poo(AGoo) 



2tt 



^p 00 (A)(G 00 ) (12) 



and the right hand side similarly. We checked the validity 
of this assumption by numerical evaluation of these two 
quantities for increasing values of n, finding less than 1% 
discrepancy for the largest clusters (n = 10 5 ). Although 
the stronger assumption, (AGoc) ~ (A) (Goo), is not valid, 




FIG. 3: The angular profile of the Jacobian factor, A n (8), 
which shows how the size of bump pre-images varies on the 
unit circle, at three stages of growth in simulations: n = 
10 3 (dashed), n = 10 4 (dash-dot) and n = 10 5 (solid). The 
exact mean-field approximation, A c (6), (thick gray) obtain 
from Eq. (7) shows a clear difference with A„(8). 



particularly near the bump center w = e l8 , the correla- 
tion seems to be canceled out in the integration along the 
angle. 

With these assumptions, we arrive at a nonlinear 
integro-differential equation for (Goo(w)) , the limiting av- 
erage cluster shape: 



(GqqH) _ Cd6 Poo (9) (|G^(e* e )|- 2 ) H e (w) 
(Gbo («?))' J^d6 Poo (6){\GU^)\- 2 ) 

7T 

d9 Poo {9) h^Heiw) 



(13) 
(14) 



where we introduce a conditional probability density, 

(\G' n (e ie )\- 2 ) 



A„(0) 



f^d0> Poo (e>)(\G> n (e ie ')\- 2 ) 



(15) 



proportional to the average local size of a bump's pre- 
image (Jacobian factor), (A n ), at angle 9. In deriving Eq. 
(13) we assume that (A^IG'J- 2 ) ~ (A 1 (t n )) (|G'J- 2 ) 
as n — > oo [16]. The stochastic nature of the aggre- 
gates is manifested through the two different averages 
in Eqs. (14)-(15). 

To check the validity of Eq. (14), we obtain from 
simulations and solve for {Goq{w)) . As shown in Fig. 3, 
the measured curves for A n (9) for n = 10 3 , 10 and 10 5 
are nearly identical, so we conclude that A 10 s (0) is a 
good approximation of A OQ (9). Now we solve Eq. (14) 
by expanding (Goo (id)) by a Laurent series and finding 
recurrence relations for the coefficients, which involve in- 
tegrals of A oc (9). We calculate 200 first coefficients and 
reconstruct {G OQ (w)) . The image of the unit circle under 
(Goo(w)), shown in Fig. 2a (thick gray), is in excellent 
agreement with the converging pattern of (G n (w)}. 

The surprising difference between the convergence 
rates of the average Jacobian A n (9) and the average map 
itself (G n (w)) is intimately related to the multifractal 
nature of the distribution of the stretching factor (har- 
monic measure) \G' n (e l& )\. Since this factor is very large 
around the cusp at 9 = 0, which is dominant during the 
growth process, fluctuations at the cusp do not contribute 
to negative moments of the distribution of \G' n (e )| and 
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thus negative moments converge much faster than posi- 
tive ones, and faster than the average map itself. Since 
A n (9) comes from averaging \G' n (e te )\~ 2 , this observation 
explains its fast convergence. This argument illustrates 
how the two averages interact in Eqs. (14)-(15) and sug- 
gests that the faster convergence of A n (9) dominates the 
morphology. 

With the validity of Eq. (14) established, we may con- 
sider its mean-field version, where the ensemble average 
is replaced by a single conformal map, given by: 



G c (w) 
A c (9) 



\ d9 Poo (9)A c (9)Hg(w), 
Jo 



J^d9 Poo (9)\G^)\- 



(16) 
(17) 



Not surprisingly, the similarity solution for continuous 
growth, Eq. (7), is an exact solution of Eqs. (16)-(17). 
In fact, it is possible to derive Eqs. (16)-(17) from a 
different representation of Eq. (6), which has been done 
for the case of DLA, Poo (9) = 1/2tt, in Rcf. [19]. Else- 
where [16], we obtain an analytical form for A c (9) for 
ADLA, which is plotted in Fig. 3 (thick gray). A small, 
but significant, difference between A c (9) and Aoo(0) is 
apparent, especially at 9 — 7r/4 and 77r/4. 

The solution in Eq. (7) can be interpreted as a self- 
consistent mean-field approximation for the average con- 
formal map, (Goo(w)) . However, fluctuations in the en- 
semble manifest themselves through the different aver- 
ages in Eqs. (14)-(15). As long as (|G-^(m;)| _2 ) is dif- 
ferent from |(GooM)T 2 > A oc (6>) ± A c (9), and thus the 
deviation of (G 00 (w)) from G c (w) is inevitable. 

We believe that the assumptions leading to Eq. (14) 
are quite general, and not specific to ADLA, so the con- 
tinuous dynamics should be a mean field theory (in this 
sense) for any stochastic aggregation, driven by confor- 
mally invariant transport processes, Eq. (1). We con- 



clude, therefore, that the solution to the continuous dy- 
namics, although similar, is not identical to the ensemble- 
averaged cluster shape. An exceptional case is DLA in 
radial geometry, where isotropy implies the trivial solu- 
tion, (G 00 (w)) — w and A c (9) = 1. Clearly, G c (w) = w 
and A c {9) = 1 solves Eq. (14) with Poc (9) = l/2vr. 

We expect, however, that this identity between the 
mean-field approximation and the average shape of 
stochastic clusters will be removed with any symmetry 
breaking, either in the model equations (such as ADLA) 
or in the BCs (e.g. DLA in a channel). This result is 
consistent with recent simulations of DLA in a chan- 
nel geometry [18], which show that the average cluster 
shape, (G n (w)), is similar, but not identical, to any of 
the Saffman- Taylor "fingers" , which solve the continu- 
ous dynamics. We expect that an analogous equation to 
Eq. (14), relating (G^w)) and (\G^{e w )\- 2 ) , will hold 
in a channel geometry, and Saffman- Taylor fingers should 
be exact solutions to the mean-field approximation of 
that equation. 

We conclude by emphasizing that, although Eq. (14) is 
a necessary condition for the average shape of transport- 
limited aggregates in the class, Eq. (1), it does not pro- 
vide a basis for complete statistical theory. Such a the- 
ory would likely consist of an infinite set of independent 
equations connecting a hierarchy of moments of the mul- 
tifractal distributions of maps {Goo(w)}, and derivatives 
{G4j(u>)}- Multifractality may speed up convergence, as 
for Eq. (14), or slow down convergence of other equa- 
tions in this set. The mean-field approximation, Eq. (16), 
which corresponds to the continuous growth process, can 
be considered as leading a hierarchy of closure approxi- 
mations to this set. 
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